Evolution and developmental expression of the sodium-iodide symporter (NIS, slc5a5) gene family: Implications for perchlorate toxicology.

The vertebrate sodium-iodide symporter (NIS or SLC5A5) transports iodide into the thyroid follicular cells that synthesize thyroid hormone. The SLC5A protein family includes transporters of vitamins, minerals, and nutrients. Disruption of SLC5A5 function by perchlorate, a pervasive environmental contaminant, leads to human pathologies, especially hypothyroidism. Perchlorate also disrupts the sexual development of model animals, including threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio), but the mechanism of action is unknown. To test the hypothesis that SLC5A5 paralogs are expressed in tissues necessary for the development of reproductive organs, and therefore are plausible candidates to mediate the effects of perchlorate on sexual development, we first investigated the evolutionary history of Slc5a paralogs to better understand potential functional trajectories of the gene family. We identified two clades of slc5a paralogs with respect to an outgroup of sodium/choline cotransporters (slc5a7); these clades are the NIS clade of sodium/iodide and lactate cotransporters (slc5a5, slc5a6, slc5a8, slc5a8, and slc5a12) and the SGLT clade of sodium/glucose cotransporters (slc5a1, slc5a2, slc5a3, slc5a4, slc5a10, and slc5a11). We also characterized expression patterns of slc5a genes during development. Stickleback embryos and early larvae expressed NIS clade genes in connective tissue, cartilage, teeth, and thyroid. Stickleback males and females expressed slc5a5 and its paralogs in gonads. Single-cell transcriptomics (scRNA-seq) on zebrafish sex-genotyped gonads revealed that NIS clade-expressing cells included germ cells (slc5a5, slc5a6a, and slc5a6b) and gonadal soma cells (slc5a8l). These results are consistent with the hypothesis that perchlorate exerts its effects on sexual development by interacting with slc5a5 or its paralogs in reproductive tissues. These findings show novel expression domains of slc5 genes in stickleback and zebrafish, which suggest similar functions across vertebrates including humans, and provide candidates to mediate the effects of perchlorate on sexual development.

Developmental tuning of mineralization drives morphological diversity of gill cover bones in sculpins and their relatives.

The role of osteoblast placement in skeletal morphological variation is relatively well understood, but alternative developmental mechanisms affecting bone shape remain largely unknown. Specifically, very little attention has been paid to variation in later mineralization stages of intramembranous ossification as a driver of morphological diversity. We discover the occurrence of specific, sometimes large, regions of nonmineralized osteoid within bones that also contain mineralized tissue. We show through a variety of histological, molecular, and tomographic tests that this "extended" osteoid material is most likely nonmineralized bone matrix. This tissue type is a significant determinant of gill cover bone shape in the teleostean suborder Cottoidei. We demonstrate repeated evolution of extended osteoid in Cottoidei through ancestral state reconstruction and test for an association between extended osteoid variation and habitat differences among species. Through measurement of extended osteoid at various stages of gill cover development in species across the phylogeny, we gain insight into possible evolutionary developmental origins of the trait. We conclude that this fine-tuned developmental regulation of bone matrix mineralization reflects heterochrony at multiple biological levels and is a novel mechanism for the evolution of diversity in skeletal morphology. This research lays the groundwork for a new model in which to study bone mineralization and evolutionary developmental processes, particularly as they may relate to adaptation during a prominent evolutionary radiation of fishes.

Pharyngeal morphogenesis requires fras1-itga8-dependent epithelial-mesenchymal interaction.

Both Fras1 and Itga8 connect mesenchymal cells to epithelia by way of an extracellular 'Fraser protein complex' that functions in signaling and adhesion; these proteins are vital to the development of several vertebrate organs. We previously found that zebrafish fras1 mutants have craniofacial defects, specifically, shortened symplectic cartilages and cartilage fusions that spare joint elements. During a forward mutagenesis screen, we identified a new zebrafish mutation, b1161, that we show here disrupts itga8, as confirmed using CRISPR-generated itga8 alleles. fras1 and itga8 single mutants and double mutants have similar craniofacial phenotypes, a result expected if loss of either gene disrupts function of the Fraser protein complex. Unlike fras1 mutants or other Fraser-related mutants, itga8 mutants do not show blistered tail fins. Thus, the function of the Fraser complex differs in the craniofacial skeleton and the tail fin. Focusing on the face, we find that itga8 mutants consistently show defective outpocketing of a late-forming portion of the first pharyngeal pouch, and variably express skeletal defects, matching previously characterized fras1 mutant phenotypes. In itga8 and fras1 mutants, skeletal severity varies markedly between sides, indicating that both mutants have increased developmental instability. Whereas fras1 is expressed in epithelia, we show that itga8 is expressed complementarily in facial mesenchyme. Paired with the observed phenotypic similarity, this expression indicates that the genes function in epithelial-mesenchymal interactions. Similar interactions between Fras1 and Itga8 have previously been found in mouse kidney, where these genes both regulate Nephronectin (Npnt) protein abundance. We find that zebrafish facial tissues express both npnt and the Fraser gene fibrillin2b (fbn2b), but their transcript levels do not depend on fras1 or itga8 function. Using a revertible fras1 allele, we find that the critical window for fras1 function in the craniofacial skeleton is between 1.5 and 3 days post fertilization, which coincides with the onset of fras1-dependent and itga8-dependent morphogenesis. We propose a model wherein Fras1 and Itga8 interact during late pharyngeal pouch morphogenesis to sculpt pharyngeal arches through epithelial-mesenchymal interactions, thereby stabilizing the developing craniofacial skeleton.

Retinoic acid metabolic genes, meiosis, and gonadal sex differentiation in zebrafish.

To help understand the elusive mechanisms of zebrafish sex determination, we studied the genetic machinery regulating production and breakdown of retinoic acid (RA) during the onset of meiosis in gonadogenesis. Results uncovered unexpected mechanistic differences between zebrafish and mammals. Conserved synteny and expression analyses revealed that cyp26a1 in zebrafish and its paralog Cyp26b1 in tetrapods independently became the primary genes encoding enzymes available for gonadal RA-degradation, showing lineage-specific subfunctionalization of vertebrate genome duplication (VGD) paralogs. Experiments showed that zebrafish express aldh1a2, which encodes an RA-synthesizing enzyme, in the gonad rather than in the mesonephros as in mouse. Germ cells in bipotential gonads of all zebrafish analyzed were labeled by the early meiotic marker sycp3, suggesting that in zebrafish, the onset of meiosis is not sexually dimorphic as it is in mouse and is independent of Stra8, which is required in mouse but was lost in teleosts. Analysis of dead-end knockdown zebrafish depleted of germ cells revealed the germ cell-independent onset and maintenance of gonadal aldh1a2 and cyp26a1 expression. After meiosis initiated, somatic cell expression of cyp26a1 became sexually dimorphic: up-regulated in testes but not ovaries. Meiotic germ cells expressing the synaptonemal complex gene sycp3 occupied islands of somatic cells that lacked cyp26a1 expression, as predicted by the hypothesis that Cyp26a1 acts as a meiosis-inhibiting factor. Consistent with this hypothesis, females up-regulated cyp26a1 in oocytes that entered prophase-I meiotic arrest, and down-regulated cyp26a1 in oocytes resuming meiosis. Co-expression of cyp26a1 and the pluripotent germ cell stem cell marker pou5f1(oct4) in meiotically arrested oocytes was consistent with roles in mouse to promote germ cell survival and to prevent apoptosis, mechanisms that are central for tipping the sexual fate of gonads towards the female pathway in zebrafish.

fras1 shapes endodermal pouch 1 and stabilizes zebrafish pharyngeal skeletal development.

Lesions in the epithelially expressed human gene FRAS1 cause Fraser syndrome, a complex disease with variable symptoms, including facial deformities and conductive hearing loss. The developmental basis of facial defects in Fraser syndrome has not been elucidated. Here we show that zebrafish fras1 mutants exhibit defects in facial epithelia and facial skeleton. Specifically, fras1 mutants fail to generate a late-forming portion of pharyngeal pouch 1 (termed late-p1) and skeletal elements adjacent to late-p1 are disrupted. Transplantation studies indicate that fras1 acts in endoderm to ensure normal morphology of both skeleton and endoderm, consistent with well-established epithelial expression of fras1. Late-p1 formation is concurrent with facial skeletal morphogenesis, and some skeletal defects in fras1 mutants arise during late-p1 morphogenesis, indicating a temporal connection between late-p1 and skeletal morphogenesis. Furthermore, fras1 mutants often show prominent second arch skeletal fusions through space occupied by late-p1 in wild type. Whereas every fras1 mutant shows defects in late-p1 formation, skeletal defects are less penetrant and often vary in severity, even between the left and right sides of the same individual. We interpret the fluctuating asymmetry in fras1 mutant skeleton and the changes in fras1 mutant skeletal defects through time as indicators that skeletal formation is destabilized. We propose a model wherein fras1 prompts late-p1 formation and thereby stabilizes skeletal formation during zebrafish facial development. Similar mechanisms of stochastic developmental instability might also account for the high phenotypic variation observed in human FRAS1 patients.

Roles of brca2 (fancd1) in oocyte nuclear architecture, gametogenesis, gonad tumors, and genome stability in zebrafish.

Mild mutations in BRCA2 (FANCD1) cause Fanconi anemia (FA) when homozygous, while severe mutations cause common cancers including breast, ovarian, and prostate cancers when heterozygous. Here we report a zebrafish brca2 insertional mutant that shares phenotypes with human patients and identifies a novel brca2 function in oogenesis. Experiments showed that mutant embryos and mutant cells in culture experienced genome instability, as do cells in FA patients. In wild-type zebrafish, meiotic cells expressed brca2; and, unexpectedly, transcripts in oocytes localized asymmetrically to the animal pole. In juvenile brca2 mutants, oocytes failed to progress through meiosis, leading to female-to-male sex reversal. Adult mutants became sterile males due to the meiotic arrest of spermatocytes, which then died by apoptosis, followed by neoplastic proliferation of gonad somatic cells that was similar to neoplasia observed in ageing dead end (dnd)-knockdown males, which lack germ cells. The construction of animals doubly mutant for brca2 and the apoptotic gene tp53 (p53) rescued brca2-dependent sex reversal. Double mutants developed oocytes and became sterile females that produced only aberrant embryos and showed elevated risk for invasive ovarian tumors. Oocytes in double-mutant females showed normal localization of brca2 and pou5f1 transcripts to the animal pole and vasa transcripts to the vegetal pole, but had a polarized rather than symmetrical nucleus with the distribution of nucleoli and chromosomes to opposite nuclear poles; this result revealed a novel role for Brca2 in establishing or maintaining oocyte nuclear architecture. Mutating tp53 did not rescue the infertility phenotype in brca2 mutant males, suggesting that brca2 plays an essential role in zebrafish spermatogenesis. Overall, this work verified zebrafish as a model for the role of Brca2 in human disease and uncovered a novel function of Brca2 in vertebrate oocyte nuclear architecture.

Sex reversal in zebrafish fancl mutants is caused by Tp53-mediated germ cell apoptosis.

The molecular genetic mechanisms of sex determination are not known for most vertebrates, including zebrafish. We identified a mutation in the zebrafish fancl gene that causes homozygous mutants to develop as fertile males due to female-to-male sex reversal. Fancl is a member of the Fanconi Anemia/BRCA DNA repair pathway. Experiments showed that zebrafish fancl was expressed in developing germ cells in bipotential gonads at the critical time of sexual fate determination. Caspase-3 immunoassays revealed increased germ cell apoptosis in fancl mutants that compromised oocyte survival. In the absence of oocytes surviving through meiosis, somatic cells of mutant gonads did not maintain expression of the ovary gene cyp19a1a and did not down-regulate expression of the early testis gene amh; consequently, gonads masculinized and became testes. Remarkably, results showed that the introduction of a tp53 (p53) mutation into fancl mutants rescued the sex-reversal phenotype by reducing germ cell apoptosis and, thus, allowed fancl mutants to become fertile females. Our results show that Fancl function is not essential for spermatogonia and oogonia to become sperm or mature oocytes, but instead suggest that Fancl function is involved in the survival of developing oocytes through meiosis. This work reveals that Tp53-mediated germ cell apoptosis induces sex reversal after the mutation of a DNA-repair pathway gene by compromising the survival of oocytes and suggests the existence of an oocyte-derived signal that biases gonad fate towards the female developmental pathway and thereby controls zebrafish sex determination.

Expression profiling of zebrafish sox9 mutants reveals that Sox9 is required for retinal differentiation.

The transcription factor gene Sox9 plays various roles in development, including differentiation of the skeleton, gonads, glia, and heart. Other functions of Sox9 remain enigmatic. Because Sox9 protein regulates expression of target genes, the identification of Sox9 targets should facilitate an understanding of the mechanisms of Sox9 action. To help identify Sox9 targets, we used microarray expression profiling to compare wild-type embryos to mutant embryos lacking activity for both sox9a and sox9b, the zebrafish co-orthologs of Sox9. Candidate genes were further evaluated by whole-mount in situ hybridization in wild-type and sox9 single and double mutant embryos. Results identified genes expressed in cartilage (col2a1a and col11a2), retina (calb2a, calb2b, crx, neurod, rs1, sox4a and vsx1) and pectoral fin bud (klf2b and EST AI722369) as candidate targets for Sox9. Cartilage is a well-characterized Sox9 target, which validates this strategy, whereas retina represents a novel Sox9 function. Analysis of mutant phenotypes confirmed that Sox9 helps regulate the number of Müller glia and photoreceptor cells and helps organize the neural retina. These roles in eye development were previously unrecognized and reinforce the multiple functions that Sox9 plays in vertebrate development.

The Fanconi anemia/BRCA gene network in zebrafish: embryonic expression and comparative genomics.

Fanconi anemia (FA) is a genetic disease resulting in bone marrow failure, high cancer risks, and infertility, and developmental anomalies including microphthalmia, microcephaly, hypoplastic radius and thumb. Here we present cDNA sequences, genetic mapping, and genomic analyses for the four previously undescribed zebrafish FA genes (fanci, fancj, fancm, and fancn), and show that they reverted to single copy after the teleost genome duplication. We tested the hypothesis that FA genes are expressed during embryonic development in tissues that are disrupted in human patients by investigating fanc gene expression patterns. We found fanc gene maternal message, which can provide Fanc proteins to repair DNA damage encountered in rapid cleavage divisions. Zygotic expression was broad but especially strong in eyes, central nervous system and hematopoietic tissues. In the pectoral fin bud at hatching, fanc genes were expressed specifically in the apical ectodermal ridge, a signaling center for fin/limb development that may be relevant to the radius/thumb anomaly of FA patients. Hatching embryos expressed fanc genes strongly in the oral epithelium, a site of squamous cell carcinomas in FA patients. Larval and adult zebrafish expressed fanc genes in proliferative regions of the brain, which may be related to microcephaly in FA. Mature ovaries and testes expressed fanc genes in specific stages of oocyte and spermatocyte development, which may be related to DNA repair during homologous recombination in meiosis and to infertility in human patients. The intestine strongly expressed some fanc genes specifically in proliferative zones. Our results show that zebrafish has a complete complement of fanc genes in single copy and that these genes are expressed in zebrafish embryos and adults in proliferative tissues that are often affected in FA patients. These results support the notion that zebrafish offers an attractive experimental system to help unravel mechanisms relevant not only to FA, but also to breast cancer, given the involvement of fancj (brip1), fancn (palb2) and fancd1 (brca2) in both conditions.

Characterization and expression pattern of zebrafish Anti-Müllerian hormone (Amh) relative to sox9a, sox9b, and cyp19a1a, during gonad development.

The role of Anti-Müllerian hormone (Amh) during gonad development has been studied extensively in mammals, but is less well understood in other vertebrates. In male mammalian embryos, Sox9 activates expression of Amh, which initiates the regression of the Mullerian ducts and inhibits the expression of aromatase (Cyp19a1), the enzyme that converts androgens to estrogens. To better understand shared features of vertebrate gonadogenesis, we cloned amh cDNA from zebrafish, characterized its genomic structure, mapped it, analyzed conserved syntenies, studied its expression pattern in embryos, larvae, juveniles, and adults, and compared it to the expression patterns of sox9a, sox9b and cyp19a1a. We found that the onset of amh expression occurred while gonads were still undifferentiated and sox9a and cyp19a1a were already expressed. In differentiated gonads of juveniles, amh showed a sexually dimorphic expression pattern. In 31 days post-fertilization juveniles, testes expressed amh and sox9a, but not cyp19a1a, while ovaries expressed cyp19a1a and sox9b, but not amh. In adult testes, amh and sox9a were expressed in presumptive Sertoli cells. In adult ovaries, amh and cyp19a1a were expressed in granulosa cells surrounding the oocytes, and sox9b was expressed in a complementary fashion in the ooplasm of oocytes. The observed expression patterns of amh, sox9a, sox9b, and cyp19a1a in zebrafish correspond to the patterns expected if their regulatory interactions have been conserved with mammals. The finding that zebrafish sox9b and sox8 were not co-expressed with amh in oocytes excludes the possibility that amh expression in zebrafish granulosa cells is directly regulated by either of these two genes.